Quantum information and entropy in random matrix ensembles 



Maciej M. Duras 

Institute of Physics, Cracow University of Technology, ulica 
Podchorazych 1, PL-30084 Cracow, Poland 

Email: mduras @ riad.usk.pk.edu.pl 

" Collective Phenomena in the Low Temperature Physics of 
Glasses"; June 16, 2003 - June 20, 2003; Max Planck Institute for 
the Physics of Complex Systems, Dresden, Germany. 

1 Abstract 

The random matrix ensembles (RME), especially Gaussian random 
matrix ensembles GRME and Ginibre random matrix ensembles, are 
applied to following quantum systems: nuclear systems, molecular 
systems, and two-dimensional electron systems (Wigner- Dyson elec- 
trostatic analogy). Measures of quantum chaos and quantum inte- 
grability with respect to eigenergies of quantum systems are defined 
and calculated. 

The distribution function for the random matrix ensembles is de- 
rived from the maximum entropy principle. Information functional 
is defined as negentropy (opposite of entropy or minus entropy). 

2 Introduction 

Random Matrix Theory RMT studies quantum Hamiltonian opera- 
tors H which are random matrix variables. Their matrix elements 
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Hij are independent random scalar variables [1, 2, 3, 4, 5, 6, 7, 8]. 
There were studied among others the following Gaussian Random 
Matrix ensembles GRME: orthogonal GOE, unitary GUE, symplec- 
tic GSE, as well as circular ensembles: orthogonal COE, unitary 
CUE, and symplectic CSE. The choice of ensemble is based on quan- 
tum symmetries ascribed to the Hamiltonian H. The Hamiltonian 
H acts on quantum space V of eigenf unctions. It is assumed that V 
is iV-dimensional Hilbert space V = ¥ N , where the real, complex, or 
quaternion field F = R, C, H, corresponds to GOE, GUE, or GSE, 
respectively. If the Hamiltonian matrix H is hermitean H = H\ 
then the probability density function of H reads: 

f H (H) = C HI3 ex P {-p-yTr(H 2 )}, (1) 

Xhh = N + ^N(N - l)f3, 
J MH)dH = 1, 

N N D-l , > 

dH=YlE n dH\]\ 
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where the parameter f3 assume values f3 = 1,2,4, for GOE(AT), 
GUE(iV), GSE(AT), respectively, and Nhq is number of independent 
matrix elements of hermitean Hamiltonian H . The Hamiltonian H 
belongs to Lie group of hermitean N x iVF-matrices, and the ma- 
trix Haar's measure dH is invariant under transformations from the 
unitary group U(AT, F). The eigenenergies E i: i = 1, A/", of i7, are 
real- valued random variables Ei = Ef. It was Eugene Wigner who 
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firstly dealt with eigenenergy level repulsion phenomenon studying 
nuclear spectra [1, 2, 3]. RMT is applicable now in many branches 
of physics: nuclear physics (slow neutron resonances, highly excited 
complex nuclei), condensed phase physics (fine metallic particles, ran- 
dom Ising model [spin glasses]), quantum chaos (quantum billiards, 
quantum dots), disordered mesoscopic systems (transport phenom- 
ena), quantum chromodynamics, quantum gravity, field theory. 



3 The Ginibre ensembles 



Jean Ginibre considered another example of GRME dropping the 
assumption of hermiticity of Hamiltonians thus defining generic F- 
valued Hamiltonian K [1, 2, 9, 10]. Hence, K belong to general 
linear Lie group GL(AT, F), and the matrix Haar's measure dK is 
invariant under transformations form that group. The distribution 
of K is given by: 

fjc(K) = C KP exp Tr(rfK)], (2) 
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Cl< > 9 = { 2ir 
M K0 = N 2 f3, 

J MK)dK = 1, 

N N D-l , > 

dK = EE E dKlj\ 

i—1 j—1 7=0 

where {3 = 1,2,4, stands for real, complex, and quaternion Ginibre 
ensembles, respectively. Therefore, the eigenenergies Zi of quan- 
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turn system ascribed to Ginibre ensemble are complex-valued ran- 
dom variables. The eigenenergies Z^i = 1, ...,iV, of nonhermitean 
Hamiltonian K are not real- valued random variables Z^ ^ Z*. Jean 
Ginibre postulated the following joint probability density function of 
random vector of complex eigenvalues Zi, Z^ for N x N Hamil- 
tonian matrices K for f3 = 2 [1, 2, 9, 10]: 

...,**) = (3) 

w 1 AT TV 

= n — 77 • n \zi - Zj\ • exp(- e Fji ), 

J=17T-J ! i<j j=l 

where Z{ are complex- valued sample points (zi G C). 

We emphasize here Wigner and Dyson's electrostatic analogy A 
Coulomb gas of N unit charges moving on complex plane (Gauss's 
plane) C is considered. The vectors of positions of charges are Zi and 
potential energy of the system is: 

U(z h ...,z N ) = - E In \zi - Zj\ + k 2 |- (4) 

i<j ^ i 

If gas is in thermodynamical equilibrium at temperature T = 7^ 
(/? = = 2, is Boltzmann's constant), then probability density 
function of vectors of positions is P(z\, Eq. (3). Therefore, 

complex eigenenergies Z{ of quantum system are analogous to vectors 
of positions of charges of Coulomb gas. Moreover, complex-valued 
spacings A 1 ^ of complex eigenenergies of quantum system: 

A 1 Z i = Z in -Z i ,i = l,...,(N -1), (5) 

are analogous to vectors of relative positions of electric charges. Fi- 
nally, complex- valued second differences A 2 Zj of complex eigenener- 
gies: 

A 2 Z t = Z l+2 - 2Z l+l + Zi,i = 1, (N — 2), (6) 
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are analogous to vectors of relative positions of vectors of relative 
positions of electric charges. 

The eigenenergies Zi = Z(i) can be treated as values of function 
Z of discrete parameter i = 1, N. The " Jacobian" of Zi reads: 

dZi A 1 Z i i 
JacZ, = — 1 ~ — 1 = A Zj. 7 
az An 

We readily have, that the spacing is an discrete analog of Jacobian, 
since the indexing parameter % belongs to discrete space of indices 
i G / = {1, N}. Therefore, the first derivative with respect to i 
reduces to the first differential quotient. The Hessian is a Jacobian 
applied to Jacobian. We immediately have the formula for discrete 
"Hessian" for the eigenenergies Z{. 

HessZ, = |J c ^§ = A 2 Zj. (8) 

Thus, the second difference of Z is discrete analog of Hessian of Z. 
One emphasizes that both "Jacobian" and "Hessian" work on dis- 
crete index space / of indices i. The spacing is also a discrete analog 
of energy slope whereas the second difference corresponds to energy 
curvature with respect to external parameter A describing parametric 
"evolution" of energy levels [11, 12]. The finite differences of order 
higher than two are discrete analogs of compositions of " Jacobians" 
with "Hessians" of Z. 

The eigenenergies Ei, i £ /, of the hermitean Hamiltonian H are 
ordered increasingly real- valued random variables. They are values of 
discrete function Ei = E(i). The first difference of adjacent eigenen- 
ergies is: 

A 1 E i = E i+l -E„i = l,...,(N-l), (9) 
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are analogous to vectors of relative positions of electric charges of one- 
dimensional Coulomb gas. It is simply the spacing of two adjacent 
energies. Real- valued second differences A 2 £^ of eigenenergies: 

A 2 E X = E l+2 - 2E i+1 + E h i = 1, (N - 2), (10) 

are analogous to vectors of relative positions of vectors of relative 
positions of charges of one-dimensional Coulomb gas. The A 2 Zi have 
their real parts ReA 2 i^, and imaginary parts ImA 2 i^, as well as radii 
(moduli) \A 2 Zj\, and main arguments (angles) ArgA 2 ^. A 2 Zi are 
extensions of real- valued second differences: 

A 2 E t = E i+2 - 2E i+1 + E i: i = 1, .., (N - 2), (11) 

of adjacent ordered increasingly real- valued eigenenergies Ej of Hamil- 
tonian H defined for GOE, GUE, GSE, and Poisson ensemble PE 
(where Poisson ensemble is composed of uncorrelated randomly dis- 
tributed eigenenergies) [13, 14, 15, 16, 17]. The Jacobian and Hessian 
operators of energy function E(i) = Ei for these ensembles read: 

dEj A^ Ei .1 . 

JacEi = —t~—^ = A Ei, 12 
oi AH 

and 

^ d Ej A Ej , o . . 

Hess ^ = ~oi * = (13) 

The treatment of first and second differences of eigenenergies as dis- 
crete analogs of Jacobians and Hessians allows one to consider these 
eigenenergies as a magnitudes with statistical properties studied in 
discrete space of indices. The labelling index i of the eigenenergies 
is an additional variable of " motion" , hence the space of indices / 
augments the space of dynamics of random magnitudes. 
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4 The Maximum Entropy Principle 



In order to derive the probability distribution in matrix space we 
apply the maximum entropy principle: 

WBx{Sf,{f x ) : (1) = 1, {H x } = Up}, (14) 

which yields: 

max{^(/^) : / f x (X)dX = 1, / H x (X)f x (X)dX = Up}, (15) 

where X = H or X = K for Gaussian or Ginibre ensembles, re- 
spectively, and TLx{X) = ^Tr(X^X). The maximization of en- 
tropy Sp(fx) = S(—kB]nfx(X))fx(X)dX under two additional 
constraints of normalization of the probability density function, and 
of equality of its first momentum and intrinsic energy, is equivalent 
to the minimization of the following functional J 7 (fx) with the use 
of Lagrange multipliers cki, Pi. 

mm{T{f x )}, (16) 
Hfx) = J(k B \nf x (X))f x (X)dX + ai j f x (X)dX 
+A / H x (X)f x {X)dX. 
It follows, that the first variational derivative of J 7 (fx) must vanish: 

^ = 0, (17) 

which produces: 

k B (ln + 1) + «i + I3{H X (X) = 0, (18) 

and equivalent ly: 

= • exp [-0 • (19) 
= exp[-(ai + 1) • /c^ 1 ], fi = Pi - k^ 1 . 
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The variational principle of maximum entropy does not force ad- 
ditional condition on functional form of 7ix(X). The information 
functional Ip is the opposite of entropy 

Wx) = Sp(f x ) = I (+k B In f x (X))f x (X)dX. (20) 
Information is negentropy, and entropy is neginformation. 
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